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Using a data sample of 106 million ^(3686) events collected with the BESIII detector operated 
at the BEPCII storage ring, we study for the first time the decay Xcj —t and Xcj —t 

K~Tv^ in the El radiative transition '0(3686) —>■ 'yXcJ- The decays are dominated by the three- 
body decay Xcj —t (pK* {892)K. We measure branching fractions for this reaction via the neutral 
and charged if* (892) and find them consistent with each other within the expectation of isospin 
symmetry. In the KKn invariant mass distribution a structure near the K* {892)K mass threshold is 
observed, and the corresponding mass and width are measured to be 1412±4(stat.)±8(sys.) MeV/c^ 
and r = 84 ± 12(stat.) ± 40(sys.) MeV, respectively. The observed state favors an assignment to 
the hi (1380), considering its possible J^^ and comparing its mass, width and decay mode to those 
reported in the Particle Data Group. 


3 


PACS numbers: 13.20.Gd, 13.25.Gv, 14.40.Pq 


I. INTRODUCTION 

It is well known that the heavy-quark mass provides 
a natural boundary between the perturbative and non- 
perturbative regimes. Quarkonium systems are regarded 
as a unique laboratory to study the interplay between 
perturbative and nonperturbative effects in Quantum 
Chromodynamics (QCD). Exclusive quarkonium decays 
constitute an important basis for investigating these ef¬ 
fects. Unlike the J/’ip and '0(3686), the P-wave char- 
monia states XcJ {J = 0, I, 2) are not directly pro¬ 
duced in e+e“ collisions, thus are less well understood 
to date [l[ . Obtaining more experimental data on exclu¬ 
sive decays of these XcJ states is important for a better 
understanding of their nature and decay mechanisms, as 
well as for testing QCD based calculations. Exclusive 
charmonium decays have been investigated widely within 
QCD. The dominant dynamical mechanism is cc quark 
annihilation into the minimal number of gluons allowed 
by symmetries followed by the creation of light quark- 
antiquark pairs, which form the final state hadrons 0. 
The XcJ states are expected to annihilate into two glu¬ 
ons. Predictions by the color singlet model give smaller 
decay widths than those determined experimentally 
[^, while much better predictions can be obtained if the 
color octet state is taken into account @,[3- Since the XcJ 
states are produced copiously in the El radiative transi¬ 
tion of '0(3686), with branching fractions around 9% [l|, 
the large 0(3686) data sample taken with the Beijing 
Spectrometer (BESIII) located at the Beijing Electron- 
Positron Collider (BEPCII) provides a unique opportu¬ 
nity for detailed studies of Xcj exclusive decays. 

In the quark model, 22 ss sector resonances, collec¬ 
tively called strangeonia, are expected below 2.2 GeV/c^. 
So far only 7 states are widely accepted experimentally, 
counting the maximally mixed rj — t]' as one ss state [8| . 
The axial-vector candidate, hi(1380), is assigned as the 
ss partner of the states, hi(1170), considering its 
mass and dominant decay to the K*{892)K final state. 
Experimentally, hi (1380) has been reported in a Par¬ 
tial Wave Analysis (PWA) only by LASS @ and Crystal 
Barrel The nature of this state is still controversial 
with respect to the predictions made by considering the 
mixing between 5'17(3)-singlet and S'C/(3)-Octet mesons 
in the l^Pi and l^Pi nonets m or those made by a 
relativized quark model M- The mass determined by 
the LASS measurement is significantly smaller than the 
theory prediction. If the LASS result is confirmed, the 
hi(1380) would seem too light to be the l^Pi ss mem¬ 
ber. The Crystal Barrel results are consistent with the¬ 
ory predictions, which means that hi (1380) could be a 
convincing candidate to be the ss partner of the l^Pi 
state hi(1170). The measurement of the total width of 
the hi (1380) is thought to be complicated by the nearby 
AT* (892)A" threshold, where the mass distribution and 


effective width can not be well described with a tra¬ 
ditional Breit-Wigner form. The direct observation of 
the hi (1380) in experiments and the precise measure¬ 
ment of its resonance parameters may shed light on its 
nature and aid in identifying the ground state axial- 
vector meson nonet in the quark model. Due to con¬ 
servation of angular momentum and parity, the axial- 
vector strangeonia candidates are not produced in J/0 
(0(3686)) radiative decays, but are expected to be pro¬ 
duced through the hadronic decay of Xcj associated with 
a vector meson 0 or in J/0 (0(3686)) decays with a 
pseudoscalar meson r]' {rj). In this paper, we report 
the first measurement of the decay XcJ ipKgK^TT^ 
and XcJ (jjK+R- 'TT^ in the electric dipole (El) radia¬ 
tive transition 0(3686) 7XcJ- In the following text, 
if not specified, AT*(892)A' denotes AT*(892)°AT*’ and its 
isospin-conjugate state Ar*(892)^Ar^, while KKtt de¬ 
notes both KgK^TT^ and AT+AT’ ■7r°. The charge con¬ 
jugated channel is always implied. This analysis is based 
on a data sample of 1.06 x 10® 0(3686) events collected 
with the BESIII detector at the BEPCII. Data with an 
additional integrated luminosity of 44.5 pb~^ [I^ at a 
center-of-mass energy of ^/s = 3.65 GeV are used to 
study continuum contributions. 


II. BESIII DETECTOR 

The BESIII detector, described in detail in Ref. 
has an effective geometrical acceptance of 93% of 47r. 
It contains a small cell helium-based main drift cham¬ 
ber (MDC) which provides momentum measurements of 
charged particles; a time-of-ffight system (TOE) based on 
plastic scintillator which helps to identify charged par¬ 
ticles; an electromagnetic calorimeter (EMC) made of 
CsI(Tl) crystals used to measure the energies of photons 
and provide trigger signals; and a muon system (MUC) 
made of Resistive Plate Chambers (RPC) located in the 
iron flux return yoke of the superconducting solenoid. 
The momentum resolution of the charged particles is 
0.5% at I GeV/c in a 1 Tesla magnetic field. The energy 
loss (dE/dx) measurement provided by the MDC has a 
resolution better than 6% for electrons from Bhabha scat¬ 
tering. The time resolution of the TOE is 80 ps (110 ps) 
in the barrel (endcaps). The photon energy resolution is 
2.5% (5%) at 1 GeV in the barrel (endcaps) of the EMC. 

Monte Carlo (MC) simulated events are used to de¬ 
termine the detection efficiency, optimize the selection 
criteria., and study possible backgrounds. A GEANT4- 
based [l^ MC simulation software, which includes the 
geometric and material descriptions of the BESIII de¬ 
tector, the detector response, and digitization models 
as well as the detector running conditions and perfor¬ 
mance, is used to generate MC samples. The 0(3686) 
resonance is simulated with the kkmc [l^ generator, 
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which is an event generator based on precise predic¬ 
tions of the Electroweak Standard Model for the process 
e+e" —>•// + ^ 7 , / = ^, T, d, u, s, c, b. The beam energy 
spread and initial state radiation (ISR) are taken into ac¬ 
count in the simulation. The subsequent decay processes 
are generated with EvtGen [l^. Background studies 
are based on a sample of 10® ^(3686) inclusive decays, 
generated with the known branching fractions taken from 
the the Particle Data Group (PDG) [I|, or with lund- 
CHARM [T^ for the unknown decays. 


III. EVENT SELECTION 

Charged particles are reconstructed from hits in the 
MDC. Charged tracks are required to be within the ac¬ 
ceptance of the MDC, satisfying | cos0| < 0.93. For each 
track, the point of closest approach to the interaction 
point (IP) must be within 1 cm in the plane perpendic¬ 
ular to the beam direction and within ±10 cm along the 
beam direction. Particle identification (PID) is carried 
out by combining information from the MDC and TOE. 
PID probabilities {prob[i)) are calculated for each track 
according to different particle hypotheses i {i = t:, K and 
p). To be identified as a kaon, a track is required to have 
prob{K) > prob{'K) and prob{K) > prob{p), while pion 
candidates are required to satisfy probiir) > prob(K) and 
probiir) > prob{p). 

Photon candidates are reconstructed from isolated 
showers in the EMC. Each photon candidate is required 
to have a minimum energy of 25 MeV in the EMC barrel 
region (|cosd| < 0.8) or 50 MeV in the endcap region 
(0.86 < |cosd| < 0.92). In order to improve the recon¬ 
struction efficiency and the energy resolution, the energy 
deposited in the nearby TOE counters is included in the 
photon reconstruction. The timing information from the 
EMC is used to further suppress electronic noise and en¬ 
ergy deposition unrelated to the event of interest. 

candidates are reconstructed with all combinations 
of two oppositely charged tracks (without a requirement 
on the point of closest approach to the IP), assuming both 
tracks to be pions. A secondary vertex fit is performed 
for each combination. Each candidate is required to 
have an invariant mass that satisfies \M^+^- — | < 

10 MeV/c^ and a decay length two times larger than its 
fit error, where M^o is the nominal mass of K'g taken 
from the PDG [l|. If more than one K'g is reconstructed 
within an event, the one with the minimum — 

M^o I is selected for further analysis. 

TT^ candidates are reconstructed from pairs of pho¬ 
tons whose invariant mass satisfies [M^o — 60] < < 

[M^o + 40] MeV/, where M^o is the nominal mass of 7r° 
taken from the PDG [l|. An asymmetrical mass window 
is used for tt® reconstruction because the photon energy 
deposited in the EMC has a long tail on the low energy 
side. A kinematic fit is performed on the selected photon 
pairs by constraining their invariant mass to the tt^ mass 


(1C fit). The Xic fro™ this kinematic fit is required to 
be less than 25. To remove backgrounds in which the 7r° 
is falsely reconstructed from a high energy photon paired 
with a spurious shower, a restriction is imposed on the 
decay angle j cosddecayj < 0.95, where ddecay is the polar 
angle of each photon in the 7r° rest frame with respect 
to the 7r° direction in the '0(3686) rest frame. If more 
than one is found within an event, the one with the 
minimum \Mry^ — M^oj is selected for further analysis. 

In the selection of the decay chain 0(3686) —>■ 7Xcj, 
XcJ 0 —>■ K'^K~, a candidate event is 

required to contain a Kg candidate, exactly four addi¬ 
tional charged tracks with zero net charge, and at least 
one photon. The four additional charged tracks must 
be identified as three kaons and one pion according to 
PID information. In the selection of 0(3686) —>■ JXcJ, 
XcJ ^ cl^K+K- 'TT^, 0 —>■ K'^K , a candidate event is re¬ 
quired to have four charged tracks with zero net charge, 
one 7r° candidate, and at least one additional photon. 
The four charged tracks must be identified as two posi¬ 
tively charged and two negatively charged kaons, respec¬ 
tively. 

To further remove potential backgrounds and to im¬ 
prove the mass resolution, a four-constraint energy- 
momentum conservation kinematic fit (4C fit) is 
performed. Events in the reaction 0(3686) —>■ 
jK+K-K^K^t:^ (0(3686) ^ xK+K-K+K-tt°) are 
required to have a xic ^ ^^0 (xic < events 

with more than one photon (besides the photons from 
the 7r° decay in the XcJ 4 >K'^K~'k^ channel), the 4(7- 
fit is repeated with each photon candidate. The photon 
candidate which gives the minimum is selected to be 
the radiative photon from the 0(3686) decay. 

In the selection of 0(3686) —>• K~K~ 

events, an additional requirement, — 

> 30 MeV/c^, is imposed to suppress the 
backgrounds 0(3686) —>■ 7r°7r°J/0 and J/0 —>• 

K~^K~K~^K~ , where Mj/^ is the nominal mass of J/0 
taken from the PDG [I|. 

After the above selection criteria are applied, the de¬ 
cay 0 —> K~^K~ is reconstructed using the two oppo¬ 
sitely charged kaons whose invariant mass is closest to 
the nominal 0 mass. Fig. [T] shows the K~^K~ invariant 
mass of 0 candidates for events in which the invariant 
mass of K'^K~KgK^TY^ (K'^K~K'^K~ is between 
3.35 and 3.6 GeV/c^. The 0 signal is described by P-wave 
relativistic Breit-Wigner function, and the background is 
a 3rd order Chebyshev polynominal function. A clear 0 
signal is observed above a very low background. Signal 
events are extracted by applying a mass window require¬ 
ment, \Mx+k- — < 15 MeV/c^, which corresponds 

to three times the mass resolution. 

The invariant mass distributions of K~ 
(K+K-K+K- 'TT*^) with the 0 mass window requirement 
are shown in Fig. [2j Significant Xcj signals are observed 
with low background. The corresponding scatter plots of 
the invariant masses of K^tt^ versus K^tt^ (K^tt^ ver¬ 
sus AT^tt*^) are shown in Fig. |31 The dominant processes 
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FIG. 1. Invariant mass distribution of K'^K for (f> candi¬ 
dates (a) for the candidate events of 'yK^K~ K%K^ TT^, (b) 
for the candidate events of K~ K~. The arrows 
indicate the (j) sideband region. 


are the Xcj —>■ (j)K*{892)K three body decays. 


IV. BACKGROUND ANALYSIS 

Since the candidate events are dominated by the 
three-body decay XcJ (j)K*{892)K, the correspond¬ 
ing branching fractions are measured by imposing a 
mass window requirement on the Ktt invariant mass, 
i.e. {Mktt — Mji*{ 892 )\ < 100 MeV/c^. For conve¬ 
nience, hereafter, an event with Xcj (pK*{892)^ 
or XcJ (/)iF*(892)°iF0 decay is called a K*{892)^ 
or AT*(892)° event. Also, an event satisfying the re¬ 
quirement \Mj^o^± — Mjf.(g 92 )±| < 100 MeV/c^ or 
\Mx±TrT—Mx*{ 892 )o\ < lOOMeV/c^ is called a Ar*(892)=*= 
or AT*(892)° tagged event. The potential backgrounds for 
the decay Xcj cpK* {892)K are studied based on the in¬ 
clusive MC sample. The following background categories 
are considered: 

The first background category contains the non- 
Xcj backgrounds, which are dominated by processes 
such as ^{8686) (/.AT*(892)7?*(892) and ^{8686) 

x4>K*{892)K . MC studies show that these backgrounds 



(GeV/c^) 


FIG. 2. Invariant mass distribution of (a) K^K KgK^n^ 
and (b) K+R-K+K' '7r° with the p mass window require¬ 
ment. 


do not produce peaks in the Xcj mass region, and their 
invariant mass spectrum can be described in the fit with 
a polynomial function when extracting the Xcj signal. 

The second background category is the non-?> back¬ 
ground, which peaks in XcJ mass region. The main 
processes of this background are XcJ /o/2'^'^*(^92)A', 
which have the same final state as that of the signal, 
where /q ^2 is ^ scalar or tensor meson, e.g. /o(1710), 
decaying to the K'^K~ final state. An unbinned maxi¬ 
mum likelihood fit is performed to the K'^K~ invariant 
mass. The fit result shows that the non-cp background is 
less than 1.8% in the (p mass window region. When we 
calculate the number of XcJ events, a (p sideband as indi¬ 
cated in Fig. [T] is used to estimate the background from 
non-?) events. The number of non-?) background events 
is subtracted to obtain the yields of Xcj signals. 

The third background category is composed of non- 
Ar*(892) events. From the scatter plots of the invariant 
mass of A'gTT^ versus AT^tt^ {K^-k'^ versus (see 

Fig. O, it appears as though the dominant backgrounds 
are the reactions containing a high mass K* state, e.g. 
Xcj ^ ?)A'*° 2 (1430)7?° and XcJ ^ <^7F*± (1430)7^^. 
This background is treated similarly to the second back¬ 
ground category: these processes have the same final 
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FIG. 3. Scatter plot of the invariant mass of (a) versus 

and (b) versus K' -7T° with the 0 mass window 

requirement. The boxes represent the IF* (892) sideband re¬ 
gions described in the text. 


state as the signal and produce peaking backgrounds in 
the XcJ mass region (background level 1.6%). Like the 
signal, XcJ —t (j>K*{892)K, this background category has 
horizontal and vertical bands around 1.430 GeV/c^, and 
contaminates in the XcJ (j)K*{892)K, illustrated as 
boxes 1 and 2 in Fig.[3l The degree of contamination can 
be evaluated using the K* 2-dimensional (2D) sideband 
region with \Mktt — 14251 < 200 MeV/c^, illustrated as 
region 3 in Fig.[3l The invariant mass for events 

within the K* 2D sideband region (without the (j) mass 
window requirement) is studied. The K^K~ invariant 
mass peaks around the cj) mass region, and the non-c^ 
events in this region are neglected in the fit. 

The last background category is the QED background, 
which is not produced via the ^(3686) resonance. Back¬ 
grounds of this type are evaluated using the data col¬ 
lected at 3.65 GeV and are found to be small and dis¬ 
tributed uniformly in the Xcj mass region. In the fit, the 
contribution from this type of background is taken into 
account by the polynomial function for the first back¬ 
ground category. 


V. SIGNAL EXTRACTION 

For the XcJ 4>KgK^T:'^ decay, the isospin con¬ 
jugate channels XcJ (/)iF*(892)°iF° and XcJ 
(j)K* {892)^ are included with K^tt^ forming a 
A*(892)° or with forming a K*{892)^. The 

branching fractions of these reactions are measured sepa¬ 
rately. The invariant mass of the K ~in dif¬ 
ferent regions is shown in Fig.|4]for (a) K*{892)^ tagged 
events, (b) iF*(892)° tagged events, (c) K*{892)^ events 
within the </> sideband region, (d) IF*(892)° events within 
the (j) sideband region, (e) events in the K* 2D-sideband 
region. 

Unbinned maximum likelihood fits are performed to 
the invariant mass spectra of K'^K ~to extract 
the yields of Xcj signal in different regions. In the fits, the 
three Xcj signals are described by the corresponding MG 
shapes convoluted with Gaussian functions which repre¬ 
sent the difference in resolution between the data and 
MG. The widths of the Gaussian functions are float. The 
background is represented by a second order Ghebyshev 
polynomial function. The results of the fits are shown 
in Fig. |4l and the signal yields are summarized in Ta- 
blelH where and are the number of 

iF*(892)± and K*{892)° tagged events, 
andiV^"':!deba„d(^*°) are those in the (j) sideband regions, 

and sideband ^-D sideband region. 

For the XcJ K~'k^ decay mode, only the 

charged IF* (892)^ state is included. In this analysis, the 
charge conjugate modes are not separated, and the joint 
branching fraction of Xcj —>■ (t>K*(892 )^is measured. 
Figure [5] shows the distributions of the K'^ K~ K'^ K~ 
invariant mass for the signal, (j) sideband, and K* 2-D 
sideband regions. 

The same fits are performed to the individual distribu¬ 
tions. The fit curves are shown in Fig. [5l and the result¬ 
ing Xcj signal yields are summarized in Table im where 
is the number of events in the signal region, 

and iV^'ldeband(^*'^) and wideband are those in the 
4> and K* 2-D sideband region. Here, e is the detection 
efficiency. 


VI. BRANCHING FRACTIONS 

The branching fractions are calculated using the ex¬ 
tracted signal yields that have been corrected according 
to the detector efficiency and decay branching fractions 
of intermediate states. 

For the XcJ decay mode, a clear in¬ 

tersection is observed on the scatter plot of the invari¬ 
ant mass of KgTT^ versus K^tt^ (Fig. E] (a)). This in¬ 
dicates that some IF*(892)=*= events are included in the 
IF* (892)° tagged events, and vice versa. Without con¬ 
sideration of the interference between two isospin conju¬ 
gated iF*(892)s in Xcj —>■ charge con¬ 

jugated iF*(892)s in Xcj —>■ (pK~^K~TT^, the relationship 
































7 



FIG. 4. The invariant mass of K~ distributions: (a) for if* (892)^ tagged events; (b) for K*{892)° tagged events; 
(c) for K*{892)^ tagged events in 0 sideband region; (d) for IF*(892)° tagged events in 0 sideband region; (e) for the events 
in the K* 2-D sideband region. 


TABLE I. The Xcj yields in different regions from fitting the K KgK^Tv^ final state. The uncertainties shown are statistical 
only. 




A'rs1deb.nd(^‘^) 


yybkg 

d) sidebandv'^ / 

^ K* sideband 

XcO 

317 ±21 

11 ± 5 

349 ± 21 

12 ±3 

14 ±4 

Xcl 

329 ± 19 

33 ±6 

310 ± 18 

22 ±5 

35 ±6 

Xc2 

443 ± 22 

21 ± 5 

428 ± 21 

17 ±4 

23 ±5 


























































TABLE 11. The XcJ yields in different regions and the detection efficiency from MC samples with K K -7r° final state. 
The uncertainties are statistical only. The statistical errors on the efficiencies are neglected because of the large size of the MC 
samples used. 




/Vrs1deband(/^**) 

jX>^kg 

K* sideband 

e 

XcO 

231 ± 17 

9±3 

9±3 

6.93% 

Xcl 

229 ± 16 

11 ±3 

18 ±4 

8.23% 

Xc2 

332 ± 19 

14 ±4 

14 ±4 

7.56% 





FIG. 5. The invariant mass distributions of K~K~ 
for the candidate events (a) within the signal region, (b) 
within the (j) sideband region, (c) within the K* 2-D sideband 
region. 


between the observed signal yields and the corresponding 
branching fractions is given by the following equations: 

TV®'® = N^(^^QgQ'jBrc€cc + n^nc^ (1) 

= N^(^2QgQ'^Brcecn + N^ip{368e)Brn£nn, ( 2 ) 


where Ai^( 3686 ) is the number of ^'(3686) events, 

is the number of iC*(892)=*=/° tagged candidates, and 
Bvc/n is the product branching fraction i?r('0(3686) —>■ 
7Xcj) X Brixcj ^ ^K*{892)^K^/^K*i892fK0) x 
Br{4> K+K-\xBr{K*{892)^ /K*{892f 

X Br{K0/K° iCg) x Br{K^ tt+tt-) 
where Br{tjj{3686) —>■ XXcj), Br{(j> —>• K^K~), 

Br{K*{892)^ if°7r±), Br{K^ j iCg), 

Br{K*(892)^ —>■ and Br{Kg —^ tt+tt”) are 

taken from the PDG [l|, and the branching fractions of 
Xcj (j)K*(892 )^/ (j)K*{892fK° are to be mea¬ 
sured in this analysis. The factor tcc/nn is the detec¬ 
tion efficiency for K*{892)^^^ events to be identified 
as AT* (892)+/° tagged events, and tcn/nc is the detec¬ 
tion efficiency for iC*(892)+/° events to be identified as 
iC*(892)+/° tagged events. Then the branching fractions 
are calculated as : 


Brc 


Brn 


fVf®e nn nc 

/Vj/’(3686)(e cc^nn ^nc^cn) 

/Ve"®ecn - iV^gecc 

/V^(3686)(e cn^nc ^nn^cc) 


(3) 

(4) 


The corresponding background subtracted signal yields 
are calculated as: 


K/l = 


- U 

- Ik* 


<"!deba„d(^^*^/°) (5) 


sideband (p sideband 
sideband 'X* sideband' 


Here, sideband and Jk* sideband are normalization fac¬ 
tors; based on the K'^K~ invariant mass distribution in 
Fig. [lia), sideband is taken as 0.37, while fx* sideband 
is taken to be 1/4. This takes into account the area 
of the sideband region (box 3 in Fig. [3|) relative to 
that of polluted signal regions (box 1 or 2 in Fig. [3]) 
as well as the existence of both isospin conjugate chan¬ 
nels in the sideband regions. The detection efficien¬ 
cies are evaluated with exclusive signal MC samples. 
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which contain events in the reactions '0(3686) —> IXcJ-, 
Xcj —t (j)K*(892)^ / (j)K*(892)^ with subsequent 
decay 0 ^ K+R-, K*{892)^ K°tt^ or K*{892)^ 
K^/K° K^, K% TT+TT-. The 0(3686) 

IXcJ decay is generated with a 1 + A cos^ 9 distribution, 
where 9 is the angle between the direction of the radia¬ 
tive photon and that of the positron beam, and A = 
1, -1/3, 1/13 for J = 0, 1, 2, assuming pure El tran¬ 
sitions. The XcJ decays are generated with a flat an¬ 
gular distribution. Due to the existence of resonances 
in the K*{892)'^K^ invariant mass, the detection ef¬ 
ficiency is evaluated by weighting the signal MC sam¬ 
ples by the invariant mass. This accounts 

for differences in the distributions between the data and 
the MC simulation. The detection efficiencies are listed 
in Table HI] and the corresponding measured branching 
fractions Br{xcj ^ (j)K*{892)^ / 4)K*{892f K°) are 
given in Table IVTIl 

TABLE III. Detection efficiencies from MC samples. The er¬ 
rors of efficiencies are neglected because of large MC samples. 



^CC 

^cn 

^nc 

^nn 

XcO 

10.47% 

4.08% 

4.29% 

10.80% 

Xcl 

11.43% 

4.34% 

4.58% 

11.96% 

Xc2 

11.42% 

4.75% 

4.49% 

11.48% 


For the XcJ 4 >K'^K~'k^ decay mode, only the 
charged /C*(892)=*= intermediate state is included, and the 
joint branching fractions of the charged conjugate modes 
are calculated by 

Brixcj ^ ^K* (892)^K^) = 

NZ - 0.25 0.5 

- U - AT -■ 

e- Br ■ iV^(3686) 

where Br is the product branching fraction of the other 
processes in the cascade decay including 0(3686) —>■ 
IXcJ, 0 —t K^K~, K*{892)^ — 7 ^ and 7 r° — 77 , 

and all the individual branching fractions are taken from 
PDG. A factor of 0.25 is determined by the K^K~ in¬ 
variant mass distribution in Fig.jTIb). The K* 2-D back¬ 
ground is taken with a weight of 0.5 due to the double 
area of box 3 in Fig. [3] (b), compared to box 1 or 2. 
As described above, the detection efficiencies are eval¬ 
uated using the exclusive MC samples weighted by the 
AT+iL“ 7 r° invariant mass. 


VII. SYSTEMATIC UNCERTAINTIES 

Several sources of systematic uncertainties are con¬ 
sidered in the measurement of the branching fractions. 
These include differences between the data and the MC 
simulation for track reconstruction, PID, photon detec¬ 
tion, kinematic fitting, selection, K'g reconstruction. 


the mass window requirement, the fitting process, back¬ 
ground estimation, MC modeling, the branching frac¬ 
tions of intermediate state decays, and the luminosity 
measurement. 

(a) Tracking and PID efficiency. The tracking effi¬ 

ciencies for and tt^ as functions of transverse mo¬ 
mentum have been investigated with the control sam¬ 
ples J/0 —>■ , Kg —>■ TT+TT” [Hi and 0(3686) —>■ 

TT+TT” J/0 [ 2 ^. The uncertainty of the tracking efficiency 
is 1% for each pion and 1.58% for each kaon. These un¬ 
certainties are obtained taking the transverse momentum 
distributions found in data into account. 

The uncertainty due to PID has been studied via the 
same control samples and is estimated to be 2 % for 
each charged pion and 2.23% for each charged kaon by 
weighting according to the transverse momentum distri¬ 
butions. 

(b) Photon detection efficiency. The uncertainty due 
to photon detection efficiency is 1% per photon [ 2 ^. This 
is determined from studies of photon detection efficiencies 
with a large and high purity control sample of J/0 —>• 

where p^ —>■ tt+tt” and tt^ —> 77 . 

(c) 4,C-kinematic fit. In the analysis, the track helix 
parameters (0o, n, tan A) are corrected for the MC sam¬ 
ple in order to reduce the difference of the 4C kinematic 
fit X 4 Q between the data and the MC sample. Here, 0o 
is the azimuthal angle that specifies the pivot with re¬ 
spect to the helix center, k is the reciprocal of the trans¬ 
verse momentum and tan A is the slope of the track. The 
correction factors are obtained from a clean sample of 
J/0 —>■ 0/o(98O), 0 —AT+AT” and /o(980) —?► tt+tt”. 
An alternative detection efficiency is evaluated with the 
same MC samples, but without helix parameters correc¬ 
tions. The difference in the efficiencies is taken as the 
uncertainty of the 4C kinematic fit psj . 

(d) TT^ selection. The uncertainty due to the tt *^ selec¬ 
tion is determined from a high purity control sample of 
J/0 —?► Tr+Tr”Tr°. The difference in the Tr° selection effi¬ 
ciency between the data and the MC simulation, 1.0%, 
is taken as the uncertainty for the Tr° selection (^ . 

(e) Kg reconstruction. The uncertainty for the A'g 
reconstruction efficiency is studied with a control sample 
of J/0 —>■ Ar*(892)+Ar+. A conservative value of 3.5% is 
taken as the systematic uncertainty for K'g reconstruc¬ 
tion 1^. This uncertainty is dominated by two sources: 
one is the tracking efficiency of the two pions from the 
ATg decay, and the other is the secondary vertex fit for 
the two pions and the related selection criteria. 

/// 0 and AT* (892) mass window requirement. The 
uncertainty from the 0 mass window requirement is es¬ 
timated by changing the 0 signal and sideband windows 
by one time the mass resolution of 0, 5MeV, which is 
determined by a fit to data in Fig. [1] . The difference in 
the branching fractions is taken as the uncertainty. For 
the uncertainty related to the AT*(892) mass window re¬ 
quirement, the value for the AT* (892) width used in the 
MC simulation, P = 50.8 ± 0.18 MeV [if, is changed by 
twice the uncertainty quoted by PDC, and the difference 
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in the detection efficiency is taken as the systematic un¬ 
certainty. 

(g) Fitting process. To estimate the uncertainties from 
the fitting process, three effects are considered. (1) Xcj 
signal lineshape. The Xcj signal lineshapes are described 
with MC simulated lineshapes convoluted with Gaussian 
functions. An alternative fit with Breit-Wigner functions 
convoluted with Gaussian functions for the XcJ signals is 
performed, where the Gaussian functions represent the 
mass resolution. The difference in the production yield 
returned by the fits is considered as the systematic un¬ 
certainty related to the signal lineshape. (2) background 
lineshape. The non-peaking background is described 
with a second order Chebyshev polynomial. Alterna¬ 
tive fits with different order Ghebyshev polynomial func¬ 
tions are performed. The largest difference in the branch¬ 
ing fractions is taken as the systematic uncertainty. (3) 
fitting range. The invariant mass of K~ 
{K~^K~K~^K~'k^) is fitted in the region of [3.35, 3.6] 
GeV/c^. Alternative fits with different ranges, [3.3, 3.6] 
GeV/c^ or [3.3, 3.65] GeV/c^ are performed. The maxi¬ 
mum difference in the branching fractions are treated as 
the systematic uncertainty. 

(h) Peaking backgrounds. The peaking backgrounds 
without a (j) signal are estimated with the events in the cj) 
sideband region [1.10, 1.13] GeV/c^. The corresponding 
uncertainty is studied by changing the (j) sideband to the 
range [1.08, 1.11] GeV/c^. The uncertainty of the peak¬ 
ing background with an excited K* is estimated by the 
difference in the branching fractions with or without the 
K* background subtraction. 

(i) Weighting method in MC. To obtain the detection 
efficiency, the MC samples are weighted by the AT* (892) A' 
invariant mass to compensate for the difference between 
the data and MC simulation. To get the uncertainty from 
this weighting method, the weight in each A'* (892) AT in¬ 
variant mass bin is randomly changed around its mean 
value by a Gaussian with a standard deviation given 
by the statistical uncertainty. The same process is per¬ 
formed ten thousand times, and the standard deviation 
on the detection efficiencies is taken as the systematic 
uncertainty related with the weighting method. 

(j) AT* (892) polarization. We estimate the uncertain¬ 
ties depending on the AT* (892) polarization by comparing 
the angular distribution of final states between data and 
MC. We treat the differences in efficiencies as the uncer¬ 
tainties. 

(k) Other uncertainties. The uncertainty of the total 
number of "0(3686) decays is 0.8%. This uncertainty is 
determined from a study of inclusive 0(3686) hadronic 
decays [2^. The uncertainties due to the branching frac¬ 
tions of intermediate states are taken from the PDG [I|. 

Table EV] and El summarize the systematic uncertain¬ 
ties for the XcJ (pKgK^iT^ and XcJ K~ 

decay modes, respectively. Assuming all of the uncer¬ 
tainties are independent, the total systematic uncertain¬ 
ties are obtained by adding the individual contributions 
in quadrature. 


VIII. THE KKtv invariant MASS AND 
111 (1380) STATE 

The KKtt invariant mass distributions are studied 
in order to identify any intermediate states. Figure [S] 
shows the distributions of the KKtt mass for the candi¬ 
date events within the Xo,i ,2 signal regions of the data 
as well as the corresponding phase space MC samples 
XcJ —t 0A'*(892)A'. The Xc 0 ,i ,2 signal regions are de¬ 
fined as [3.365, 3.455] GeV/c^,’[3.490, 3.530] GeV/c^ and 
[3.540, 3.575] GeV/c^, respectively. A threshold enhance¬ 
ment, which can not be described with the phase space, 
is observed in both Xci and Xc 2 signal regions (Fig.[6](b), 
(c), (e), (f)), but is absent in the Xco signal region (Fig. [5] 
(a), (d)). 

Assume that the threshold enhancement comes from 
an intermediate state, called X. Since the X state is pro¬ 
duced in the Xci ,2 decay associated with a 0 meson, its 
C parity is must be negative due to the conservation of 
C parity. The X state is observed in the A'* (892) A final 
state, and in the Xc 2 decay process associated with a 0 
meson, which rules out an assignment of = 0+. If 
we assume that the X state is a conventional qq state, we 
can also rule out a of 0“, 2+ etc. Taking into account 
its mass region, its decay through Ar*(892)Ar^nd the ac¬ 
cessible and comparing with the PDG [II], the most 
likely assignment for X is the /ii(1380) = I"' ). In 

the following, we assume that the AT* (892)10 threshold 
enhancement is the hi (1380) intermediate state. 

Besides the AT* (892)A threshold enhancement ob¬ 
served in the Xci ,2 signal regions, a clear structure around 
1.7 GeV/c^ is observed in the XcO decay, also evident in 
the Xci ,2 decay (Fig. [5]). Gonsidering the mass region, 
and its decay through AT* (892) A, we conclude that this 
structure is the 0(1680). The C parity should be minus, 
and the 0(1680) meets this requirement. There is also a 
possible 0(1850) contribution visible in Fig. [B] 

To extract the resonance parameters of the hi(1380), 
a simultaneous fit is performed to the ATgAT^Tr^ and 
invariant mass distributions for the Xci ,2 can¬ 
didate events (Fig. [5] (b), (c), (e), (f)). Three reso¬ 
nance states, hi(1380), 0(1680) and 0(1850) as well as 
a phase space contribution from Xci ,2 —t 0Ar*(892)A' are 
included in the fit without interference. The phase space 
distribution is described using the shape of the MG sam¬ 
ples, smoothed with the KEYS method [l^. The shapes 
of 0(1680) and 0(1850) are described by relativistic Breit- 
Wigner functions with constant widths. All theparame- 
ters of the resonances are fixed to PDG values [l|. Based 
on the results of the Xcj fits, the background from non- 
Xcj contributes is about 2.7%. MC studies show that 
the shape of the KKtt invariant mass of the non-Xcj 
contribution is similar to that of the phase space. There¬ 
fore, this background is included in the phase space con¬ 
tribution in the fit to the K^K^tt^ and K^K~7t^ in¬ 
variant mass distributions. Since the /ii(1380) resonance 
is close to the AT* (892)A' mass threshold, it is parame¬ 
terized with a relativistic Breit-Wigner function with a 
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TABLE IV. Systematic uncertainties on the branching fraction for the XcJ 




final states (in %). 


Sys. err source 

XcO 

Xci 

Xc2 

Common contribution 

Total number of '1/^(3686) 


0.8 


Tracking 


5.7 


Particle ID 


8.7 


Photon selection 


1 


Kinematic Fit 


1 


reconstruction 


3.5 


(p mass window 


1 


Xcj ^ 0A*(892) 

(xcj 

(j)K*{892)'^RO) 


K* mass window 

1.0 (2.1) 

0.5 (0.6) 

1.0 (2.1) 

XcJ lineshape 

1.5 (3.4) 

2.3 (3.0) 

0.0 (3.7) 

Fit range 

2.0 (4.2) 

1.1 (1.2) 

0.7 (1.4) 

Non-peaking BG shape 

1.0 (0.0) 

0.5 (0.0) 

0.7 (0.5) 

Peaking BG without (p 

0.0 (0.5) 

0.6 (0.7) 

0.8 (0.0) 

Peaking BG with (p 

1.0 (0.8) 

3.2 (3.0) 

1.3 (1.4) 

BRs from PDG 

2.7 (2.7) 

3.3 (3.3) 

3.4 (3.4) 

MC model 

2.7 (2.6) 

2.3 (2.8) 

2.1 (2.2) 

K*(892) polarization 

6.1 (5.0) 

7.1 (6.5) 

6.6 (5.7) 

Sum 

13.6 (14.0) 

14.5 (14.4) 

13.7 (14.0) 


TABLE V. Systematic uncertainties on the branching fraction for the Xcj —t 4>K'^ K 7r° final states (in %). 


Sys. err source 

XcO 

Xci 

Xc2 

Common contribution 

1 / 1 ( 3686 ) total number 


0.8 


Tracking 


6.3 


BID 


8.9 


Photon selection 


1x3 


7r° reconstruction 


1 


Kinematic Fit 


1 


(j) mass window 

0.5 

1.2 

2.2 

K* mass window 

0.5 

2.7 

2.6 

Xcj lineshape 

8.5 

9.3 

6.2 

Fit range 

2.4 

1.1 

1.0 

Non-peaking BG 

3.3 

1.6 

0.7 

Peaking BG without (j> 

0.5 

0.6 

0.4 

Peaking BG with (j) 

1.9 

4.3 

2.0 

BRs from PDG 

2.7 

3.3 

3.4 

MG model 

2.0 

1.0 

0.9 

A* (892) polarization 

6.6 

6.2 

7.4 

Sum 

16.8 

17.3 

15.9 


mass dependent width . A detailed description of the 
hi (1380) lineshape used in the following fits is given in 
Appendix The simultaneous fit to the K^K^-k^ and 
K'^K~'k^ invariant mass distributions is performed for 
the candidate events in the Xci a-nd Xci mass regions. 
The fit yields a mass of (1412.4 ± 4.4) MeV/c^ and a 
width of (84.5 ± 12.4) MeV for the hi (1380) resonance, 
where the errors are statistical only. In the fit, we take 
the change of the detection efficiency as a function of the 
invariant mass into account. The efficiency functions are 
0.12 - 0.02 X , 0.14 - 0.03 x , 0.18 - 

0.02 X Mx+K-Tr° and 0.18 — 0.03 x Mx+k--k° for the Xci 
and Xc 2 regions, respectively. The product of the rela¬ 


tivistic Breit-Wigner functions and the efficiency func¬ 
tions are used to describe the signal. Figure [7] shows the 
sum of the invariant mass distribution of the 
and decay modes for the candidate events in 

the Xci and Xc 2 mass regions as well as the sum of the cor¬ 
responding fit curves. The goodness of fit is determined 
to be x^/n.d./ = 1.09 by projecting all candidate events 
in 45 bins. The statistical significance of the hi (1380) 
signal is measured to be greater than 10 tr by compar¬ 
ing the likelihood values of the fit with and without the 
hi (1380). The statistical significances of the ())(1680) and 
())(1850) signals evaluated by the same method are found 
to be 4.3(7 and 3.2 (t, respectively. We consider the inter- 
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FIG. 6. The KK-k invariant mass distributions, (a) in the XcO region; (b) in the Xci region; (c) 

in the Xc 2 region; (d) K~tv° in the XcO region; (e) K~tx^ in the Xci region; (f) K~tt^ in the Xc 2 region. The black 

dots represent the data, and the histogram shows the phase space MC sample. 


ferences between /ii(1380) and these resonances in the 
systematic uncertainties. 

The different sources of systematic uncertainty for the 
hi (1380) resonance parameters are considered as follows: 

(a) Parameterization of the energy dependent width. 
The energy dependent width of the hi (1380) is described 
with the truncation functions combined of exponential 
and polynomial functions. The curves of S{m) are re¬ 
fitted with the same truncation functions, but with dif¬ 


ferent truncation point. A fit was also attempted to the 
S(rn) curves using a second-order polynomial. The new 
fitted functions of S{m) are used in the simultaneous fit 
on the KKtt invariant mass. The largest differences in 
the mass and width of the hi(1380), 3.9 MeV/c^ and 3.8 
MeV, are taken as one of systematic uncertainties for the 
hi (1380) resonance parameters. 

(b) Background shape. In the fit, the background is 
described by a smoothed phase space MC sample. An al- 




























































































































































































































13 



FIG. 7. The sum of and K^K^-k° mass spectra 

in the Xci and Xc 2 mass regions. The markers with error 
bars represent the data; the dash curve the /ii(1380) signal; 
the dash-dot-dot curve the <j!)(1680) signal; and the dash-dot 
curve the (?i(1850) signal. 


ternative fit is performed using a second-order Chebychev 
polynomial function to describe the non-resonant. The 
differences in mass, 4.7 MeV/c^, and width, 5.5 MeV, are 
taken as the systematic uncertainties due to the shape of 
the background. 

(c) Fit range. An alternative fit with a different fit 
range, [1.2, 2.0] GeV/c^, is performed. The differences in 
the mass and width of the /ii(1380) are found to be 0.5 
MeV !(? and 3.5 MeV, respectively. 

(d) Efficiency parameterization. A fit with a flat 
efficiency curve is performed, and the differences of 0.3 
MeV/c^ in the mass and 1.0 MeV in the width are taken 
as the systematic uncertainties related to the efficiency 
correction. 

(e) 0(1680) resonance. In the nominal fit, the mass 
and width of the 0(1680) resonance are fixed to their 
PDG values. Different fit scenarios, (1) without the 
0(1680), (2) leaving the mass and width of the 0(1680) 
resonance free in the fit, and (3) changing the mass and 
width by one standard deviation from the PDG values are 
performed, and the largest change in the mass and width, 
1.4 MeV/c^ and 17.0 MeV, are regarded as the system¬ 
atic uncertainties related to the 0(1680) resonance. 

(f) 0(1850) resonance. The systematic uncertainty 
due to the 0(1850) resonance is evaluated in a similar way 
as for the 0(1680). The largest change in mass and width, 
1.1 MeV/c^ and 3.5 MeV, are taken as the systematic 
uncertainties related to 0(1850) resonance. 

(g) The branching fraction constraint between isospin 

conjugate decay modes. In the nominal fit, no con¬ 
straint is imposed on branching fractions for the the 
isospin conjugate modes and different final states. The 
systematic uncertainty related to the branching fraction 
constraint is estimated by requiring that the branch¬ 
ing fraction of /ii(1380) —>■ is twice that of 

/ii(1380) —>■ K^K~'k^, which is expected by isospin 
symmetry. The changes in the mass and width of the 


/ii(1380), 0.3 MeV/c^ and 3.4 MeV, are taken as the 
systematic uncertainties of this type. 

(h) Interference effects. The fit is repeated under 
the following interference scenarios, (1) the interference 
between /ii(1380) and phase space, (2) the interference 
between ft.i(1380) and 0(1680), (3) the interference be¬ 
tween ft.i(1380) and 0(1850), The largest differences to 
the nominal values of the mass and width, 4.7 MeV/c^ 
and 35.3 MeV, are taken as the systematic uncertainties 
related to interference effects. 

Table EH shows the systematic uncertainties for the 
mass and width of hi (1380). Assuming that all sources 
of systematic uncertainty are independent, the total sys¬ 
tematic uncertainty is determined to be 8 MeV jc? for the 
mass and 40 MeV for the width. 


IX. SUMMARY AND CONCLUSIONS 

Using 106 M 0(3686) events collected with the BE- 
SIII detector, we present a study of the decays XcJ 
(fKgK^TT^ and XcJ 4>K'^K~n^, via the El radia¬ 

tive transition 0(3686) —>■ xXcJ- The decays are dom¬ 
inated by the three-body reaction XcJ (j)K*{892)K. 
The branching fractions for this reaction via neutral and 
charged A* (892) are measured for the first time and are 
summarized in Tables ETIl The measured branching frac¬ 
tions for XcJ —>■ 0Ar*(892)^ArT jn the two different final 
states are found to be consistent with each other. The 
ratio of the branching fraction_for Xcj —>• 4>K*(892)^ 
to that of Xcj —>■ 0Ar*(892)°iT° is consistent within the 
expectations of isospin symmetry. 

By examining the invariant mass spectrum of KKn, 
a significant excess of events above the phase space ex¬ 
pectation is observed near the AT* (892)A mass thresh¬ 
old in the decays of Xci,2, with a significance greater 
than 10(7. The observed structure has negative C par¬ 
ity, and is expected to be the ft.i(1380) state, consider¬ 
ing its mass, width and decay through K*{892)K. A 
simultaneous fit is performed to the invariant mass dis¬ 
tributions of KKtt for the candidate events in the Xci,2 
signal regions. The mass and width of the /ii(1380) 
are determined to be 1412 ± 4(stat.) ± 8 {sys.) MeV/c^ 
and 84 ± 12{stat.) ± 40(s?/s.) MeV, respectively. This 
is the first direct observation of the /ii(1380) in its de¬ 
cay to K*{892)K. Evidence is also found for the decays 
XcJ 00(1680) and XcJ 00(1850), but with signif¬ 
icances less than 5(7. More data and advanced analysis 
techniques, e.g. PWA, may shed light on the properties 
of the structures observed in the KKtt invariant mass 
spectrum. 
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TABLE VI. Systematic uncertainties on mass and width for the hi(1380). 


Sys. err source 

mass (MeV) 

width(MeV) 

Parameterization of energy dependence width 

3.9 

3.8 

Shape of the background 

4.7 

5.5 

Fit range 

0.5 

3.5 

Efficiency curve 

0.3 

1.0 

(;i(1680) resonance 

1.4 

17.0 

0(1850) resonance 

1.1 

3.5 

Isospin constraint 

0.3 

3.4 

Interference effects 

4.7 

35.3 


Sum 


7.9 40.2 

TABLE VII. 

Branching fractions measured in 

(pKK-K final states 

Decay Modes 

(PKsK^tv^ (x10“^) 

<j>K+K-Tr'^ (xlO"^) 

0A'*(892)±A:=f 
0A*(892)°AO-bc.c. 

1.65 ± 0.21{stat.) ± 0.22{sys.) 
2.03 ± 0.21(sfat.) ± 0.28(sys.) 

1.90 ± 0.14(stat.) ± 0.32{sys.) 

0A'*(892)±A=F 
0A*(892)°A:o- bc.c. 

1.76 ± 0.21(sfat.) ± 0.26(sys.) 
1.51 ± O.wlstat.) ± 0.22(s2/s.) 

1.62 ± 0.12(stat.) ± 0.28{sys.) 

0A'*(892)±A:+ 

(j)K* (892)°RO + C.C. 

2.56 ± 0.23(stat.) ± 0.35{sys.) 
2.27 ± 9.22^stat.) ± 9.32{sys.) 

2.74 ± 0.16{stat.) ± 0.U{sys.) 
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Appendix A: Appendix: Lineshape of hi (1380) 

The hi (1380) resonance is parameterized with a rel¬ 
ativistic Breit-Wigner function with a mass dependent 


width 

TOgr^(m;mo) 

(mg — m2)2 + mQr2(m; mg) ’ 

where mo is the nominal mass of hi (1380) state, and 
F(m;mo) is the corresponding mass dependent width, 
which will be discussed in the following. 

In typical two-body decays, the width r(m, mg) varies 
with mass roughly as [2^ 

r(m, mo) ~ ro(-^)^'+\ (A2) 

90 

where I is the orbital angular momentum, q is the mo¬ 
mentum available in a normal two-body decay, and Fo 
and go are the corresponding widths and momenta in 
the nominal mass. However, in a quasi two-body decay 
such as hi(1380) — K*{892)K, the q of the decay iso¬ 
bar (AT*(892)) is no longer precisely defined. Due to the 
non-zero AT* (892) width, the A'*(892)A' threshold is not 
well defined and the q momentum available in a normal 
two-body amplitude may become un-physical {q^ < 0). 
Following Ref. [s^, if we assume the energy dependence 
of F is proportional to the integral of the decay matrix 
element over the available area of the Dalitz plot for each 
KKtt mass value, then 
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r(TO;TOo) = To 


S{m) 

S'(mo) 


Dalitz plot{m) 

7 


Dalitz pZot(mo) 


(A3) 


and for S'-wave AT* (892)AT 

-Dk* = I BWki^ ■h+9- BWki, •^'*21, ( A4) 

where BWk*. are the usual P-wave relativistic Breit- 
Wigner functions [2^ with the mass and width of 
AT* (892) fixed at their PDG values. Here U are the 
Zemach vectors describing the spin 1 content, and are 
the vector part of 

9 2 

TJl^ _ TYl,^ 

(A5) 


and Pi (i = 1, 2, 3) are the vector parts of the 4-momenta 
of TT, K and K in the KKtt rest system. The interference 
between the K*{892)K and K*{892)K is fixed by C par¬ 
ity: constructive {g = -fl) for C = -|-1, and destructive 
{g = —1) for C = —1. In our analysis, C = —1. 


To simplify the fitting process, a numerical method is 
used to extract the integral of the decay matrix element 
over the available phase space. The integral can be re¬ 
written as : 


S{m) = 


n2 

Uk 


Dalitz plot{m) 


2=1 2=1 


(A6) 


where the term (i) can be extracted using 

the MC samples of Xcj —t 4>2C (A —>■ KKn) with a large 
statistics N, and V is the phase space volume for the 
decay X —>■ KKtt. 

A set of MC samples generated with a different mass 


m and zero width for the intermediate state X is used to 
extract (0- The phase space volume V of 

a three-body decay with system mass m can be derived 
as [3l|: 


R'i{m\mi,m2,m^) = / R2{m\mi2,m^)R2{rnl2-,rni,rn2)dm\2 

J (mi+m2)^ 


r 


(A7) 


where mi2 is the invariant mass of system of particle 

I 


1 and particle 2, and R 2 is the two-body phase space 
volume, which can directly taken from PDG [l| : 


i?2(mi2;mi,m2) = J[l- (I!^l±^)2][l - (A8) 

V TO12 mi2 


The above integral process gives the same result as the 
integral over the decay area in the Dalitz decay study . 

The value of ^ the phase space volume 

V, and the corresponding product, S{m), as function of 
invariant mass of KKtt system, m, are shown in Fig. [S] 
Due to the small difference in the mass and width of neu¬ 


tral and charged A*(892) and Xci ,2 phase space volume, 
the calculations are performed for Xci ,2 and the decay 
mode of K^K^tt^ and K^K~tt^, individually. S{m) 
represents the mass dependence of the width r(m,m'^) 
for the hi(1380), which is used to describe the hi(1380) 
line-shape. 
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